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I. INTRODUCTION 

The successful explanation of the space-time evolution of matter produced in relativistic heavy-ion collisions at 
RHIC (Relativistic Heavy-Ion Collider) and the LHC (Large Hadron Collider) using relativistic viscous hydrody¬ 
namics initiated new developments in the hydrodynamics [1-24]. An example of such an activity is the burgeoning 
study of anisotropic hydrodynamics. In its original formulation, anisotropic hydrodynamics was restricted to boost- 
invariant and transversally homogeneous systems [25, 26], denoted below as the (0+l)D case. Shortly afterwards, this 
approach was generalised to the one-dimensional, non boost-invariant situation [27, 28]. The (2+l)D and (3+l)D 
cases (corresponding to a two-dimensional transverse expansion with longitudinal boost-invariance and to an arbitrary 
three-dimensional expansion without any symmetry constraints, respectively) were studied in a subsequent series of 
papers [29-32] . In parallel, the concept of anisotropic hydrodynamics was extended to describe mixtures of anisotropic 
quark and gluon fluids [33-35] . 

In Refs. [25-35] only the leading order of the hydrodynamic expansion has been included and described by the 
spheroidal Romatschke-Strickland form [36]. This form allows for the difference between the longitudinal and trans¬ 
verse pressures only (with the cylindrical symmetry in the transverse plane) and hence cannot describe more complex 
situations where, for example, the three components of pressure differ from each other. A successful method to solve 
this problem by including further terms in the hydrodynamic expansion was first presented in Ref. [37] for massless 
(conformal) systems and was recently extended to massive (non-conformal) systems in Ref. [38]. 

An alternative starting point for performing the anisotropic hydrodynamics expansion is to include the momentum 
anisotropy in a more complete way already at leading order by generalizing the Romatschke-Strickland form to allow 
the energy-momentum tensor to possess different pressures along all three spatial directions in the local rest frame 
(ellipsoidal form). Such a non-perturbative treatment for the (1+1 )D case (boost-invariant expansion with cylindrical 
spatial symmetry) was proposed for conformal systems in Ref. [39] and soon generalized to non-conformal systems in 
Ref. [40]. Recently, a more general set of equations for leading-order anisotropic hydrodynamics has been presented 
by Tinti in Ref. [41]. Differently from previous leading-order approaches, the new formulation does not use any 
simplifying symmetry assumptions such as the longitudinal boost invariance or cylindrical symmetry. 

Using this new approach it was shown that, in the close-to-equilibrium limit, the general anisotropic hydrodynamics 
equations have the same structure as the corresponding second-order viscous hydrodynamics equations [41], However, 
in the process one finds that anisotropic hydrodynamics allows for several different ways to obtain the equations 
necessary to describe the evolution of the bulk pressure of the system. In order to determine which option is the 
best for possible phenomenological applications, in this paper we test the corresponding solutions with the help of 
exact solutions of the Boltzmann equation treated in the relaxation time approximation. To perform the numerical 
comparisons, the new equations are reduced to the boost-invariant and transversally homogeneous case, where the 
solutions of the Boltzmann equation are available. We note that comparisons of the predictions of various hydro- 
dynamic frameworks with the solutions of the underlying kinetic theory has now become a quite useful method for 
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testing the accuracy of different hydrodynamic formulations [42-48]. 

In this paper, we show that the agreement of the new formulation of anisotropic hydrodynamics with the exact 
solutions is similar to that obtained using earlier, less general, formulations of anisotropic hydrodynamics. We find 
that leading-order anisotropic hydrodynamics better reproduce the exact solution for the pressure anisotropy and gives 
comparable results for the bulk pressure evolution. With regard to which method to use to describe the evolution of 
the bulk viscous pressure, we find that it is better to use the zeroth moment of the kinetic equation, while for the 
ratio of the longitudinal to transverse pressure it is best to use the trace of the second moment instead of the zeroth 
moment. 

The paper is organized as follows: In Sec. II, we present the equations of leading-order (3+l)D anisotropic hydro¬ 
dynamics derived originally in Ref. [41] — we first define the energy-momentum tensor in the Landau frame, then the 
phase space distribution is introduced and discussed, and finally, the forms of the hydrodynamic equations obtained 
from the first three moments of the Boltzmann equation are presented. In Sec. Ill, the symmetry conditions character¬ 
izing boost-invariant and transversally-homogeneous systems are implemented and special forms of various physical 
quantities in the (0+l)D case are presented. In Sec. IV, we present the leading-order anisotropic hydrodynamics 
equations reduced to the (0+l)D case. In Sec. V, we present our numerical results, first for the fully non-perturbative 
case (Sec. VA), and then for the linearized case (Sec. VB). We summarize and conclude in Sec. VI. 


II. LEADING ORDER (3+l)D ANISOTROPIC HYDRODYNAMICS 
A. Energy-momentum tensor 

In this Section we recall the general structure of the equations of (3+l)D leading-order anisotropic hydrodynamics 
as introduced recently by Tinti [41]. In this case, one adopts the Landau frame and, consequently, uses the following 
form of the energy-momentum tensor 

T^ v = Tg” + — IIA F = £U fl U v -(v eq . + n) A''" + Td 11 '. (1) 

Here the tensor Tg” describes the equilibrium part of , which has the perfect-fluid form 

Tg = £eq .U^ir - Peq.A^, (2) 

expressed by the equilibrium energy density, 

£ eq . = 4 t TNT 4 ml q [: 3K 2 {rh eq .) + m eq .A'i(m eq .)], (3) 

the equilibrium pressure V eq ., 

V eq . = 4t tN T A m | q . K 2 (m eq .), (4) 

and the flow vector tVL In Eqs. (3) and (4) N = A^of/(27r) 3 , with Vdof being the number of internal degrees of 
freedom, m eq . = m/T, where m is the particle mass and T is the temperature. We note that Eqs. (3) and (4) 
describe a classical massive gas using a Boltzmann distribution. Such a system has been recently studied also in 
Refs. [40, 49, 50]. Generalizations to the case of quantum (Bose-Einstein or Fermi-Dirac) statistics case can be found, 
e.g., in Refs. [51, 52]. 

The projection tensor A**" is defined as 

A**" = - U»U V . (5) 

The quantity n is the bulk viscous pressure, while 7is the shear tensor. For conformal systems (in particular for a 
gas of massless particles) the bulk viscous pressure vanishes. On the other hand, for non-conformal systems, the bulk 
viscous pressure is finite and there are substantial couplings between the bulk pressure and shear tensor [47, 48]. 


B. The underlying phase-space distribution function and Landau matching 


The energy-momentum tensor T is the second moment of the phase-space distribution function f(x,p), which 
appears if we analyze the first moment of the Boltzmann equation, 

rj-iflU _ 


dPp^p v f{x,p). 


( 6 ) 
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Here dP is the Lorentz-invariant momentum integration measure dP = d 3 p/E p , where E p = \/p 2 + m 2 is the particle 
energy. In the present formulation of anisotropic hydrodynamics we assume that the distribution function is of the 
form [30, 40] 


f = N exp 



( 7 ) 


Here, the tensor reads [40] 


= U»U V + ^ v - cj> A'"', (8) 

where is the spatial, traceless part of 'EJ iv , and <j> is a degree of freedom controlling the trace of 3 f “'. The expression 
(7) is a generalized Romatschke-Strickland form [36]. 

The equilibrium distribution function is recovered in the limit of vanishing anisotropy tensor and vanishing 
bulk parameter <j>. In this case, one has 


/ = fe q . = N exp 


T 


( 9 ) 


For any non-equilibrium phase-space distribution function /, one can find the corresponding thermal distribution / eq . 
by requiring that the two distributions yield the same energy density, namely 

£ = / dP(p ■ U) 2 f(x,p) = J dP(p- U) 2 f eq .(x,p) = £eq.- (10) 

Equation (10) is known as the Landau matching. It allows us to determine the effective temperature of the system 
locally for any non-equilibrium situation. We note that in local equilibrium, the momentum scale A is equivalent to 
the temperature T. 


C. First-moment equations 


Similarly to other hydrodynamic prescriptions, the first four equations of (3+l)D anisotropic hydrodynamics follow 
from energy-momentum conservation 


d tl T' lv = 0. (11) 

We call these equations the first-moment equations as they are obtained from the first moment of the Boltzmann 
equation. In the relaxation time approximation, Eq. (11) follows from the Landau matching condition (10). 

Instead of using (11) directly, we consider the projection of (11) on the four-velocity U, 

D£ e q. = - (£eq. + V eq . + n) e + 7(12) 

and the contraction of (11) with the projector A, which gives 


(feq. + V eq . + n) DU» = V"p eq . + v M n + A^d^. (13) 

In Eq. (12) the symbol D = is the convective derivative (i.e., the time derivative along the fluid flux lines), 6 

is the expansion scalar 9 = d^U^, and a^ = d l ' IJ U lJ ' 1 is the velocity stress tensor. The angular brackets denote the 
following tensor structure 

A <^> = = 1 ^ A £ A g + _ 'l A ^A af) 

The operator V /; appearing on the right-hand side of Eq. (13) is the spatial gradient V^ 1 = A^d^. The latter can be 
used in the standard decomposition of the four velocity gradient 



^iiUu — &[iif T ^dA^ T 


(15) 


where u> pl , = — V^C/^) = is the vorticity. 









4 


D. Second and zeroth moment equations 

The analysis of the second moment of the Boltzmann equation performed in [41] shows that it is convenient to use 
the tensor 

0^ = JdP(p ■ U)p a pPf. (16) 

By taking traceless part of 0^", we obtain the tensor 0^©, which satisfies the following differential equation [41] 

+ — 0^"> = -20<>© A - -0 0^© + 20jfV'>\ (17) 

7"eq. 3 

where r eq . is the relaxation time appearing in the Boltzmann kinetic equation written in the relaxation time approxi¬ 
mation [53]. We note that, in the close-to-equilibrium case, the tensor 0^© becomes proportional to the shear tensor 
7 . Besides the tensor equation (17), which describes effects related to the shear viscosity, in this paper we will 
consider three prescriptions for obtaining the additional scalar equation necessary to describe the evolution of the 
bulk parameter cj). 


1. Trace equation 

The first possibility we consider results from taking the trace of the second moment of the Boltzmann equation 
which gives [41] 


£>0 tr + |6l0tr-20 <M ^ = — 


0t° r q ' ^ ©t 


(18) 


where 


©tr = -A a/3 0^. 


(19) 


2. U projection 


The second possibility we consider is obtained by projecting the second moment of the 
times on the flow vector U. In this way one obtains [41] 


Boltzmann equation two 


m 2 (^Dn + nd'j + 


£>0 tr + -00tr-20 <M a' 


p V 


m 


• eq. 


n eQ . - n 


+ 


0 t eq - - 0 t 


( 20 ) 


where n eq , and n denote the particle number density in and out of equilibrium, respectively. Note that Eq. (20) 
reduces to Eq. (18) in the limit m —> 0. 


3. Zeroth moment equation 

The third option we consider is to use the zeroth moment of the Boltzmann equation instead of the second moment. 
In this case, one obtains 


Dn + n9 = 




( 21 ) 


It is interesting to observe that Eq. (20) is the sum of Eq. (18) and Eq. (21) multiplied first by m 2 . We note, that 
in general, any linear combination of Eqs. (18) and (21) is acceptable as a scalar equation describing the bulk effects, 
not only Eq. (20). 

Before proceeding further, we would like to discuss the independent degrees of freedom present in our approach. 
Since all physical variables are expressed by the non-equilibrium distribution function (7), they should depend on: 
the momentum scale parameter A, three independent components of the flow vector U, five independent components 
of the tensor . and the parameter q !>. These are altogether ten independent variables. Thus, in the (3+l)D case 
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we need ten dynamic equations, which can be taken as Eqs. (12), (13), (17) and one equation out of Eqs. (18), (20), 
and (21). 

We note that the effective temperature T is always calculated from the Landau matching (10). If we want to treat 
T as an independent variable, we have to consider the Landau matching as an additional constraint needed to close 
the system of equations. 


III. IMPLEMENTATION OF THE BJORKEN SCENARIO 
A. Local orthonormal basis for the Bjorken flow 


In the Bjorken flow limit, the symmetry of the evolution constrains and greatly simplifies the equations of motion. 
Because of boost invariance in the longitudinal direction and rotational invariance in the transverse plane, the flow 
four-velocity I/ M reads 


( 22 ) 

where ry is the space-time rapidity. Each scalar function depends in this case only on the longitudinal proper time 

t = \/t 2 — z 2 . (23) 

In the local rest frame there is only one privileged space direction, along the beam axis, namely 

Z* = (sinhry, 0,0, cosh V ), Z ■ Z = -1, Z-U = 0. (24) 

In order to have a complete orthonormal basis, we can choose the other two four-vectors X and Y in such a way that 
they span the transverse plane and are orthogonal to the U and Z vectors, namely 

XX = -1, Y-Y = - 1, X ■ U = 0, XZ = 0, Y-U = 0, Y ■ Z = 0. (25) 


U M = (cosily, 0, 0, sinhry), 


?/ =-ln 


t + z 
t — z 


They can be simply the directions of the x and y coordinates in the lab frame. Therefore, seen from the local rest 
frame (LRF), they read 


£7" = (1,0,0,0), X" = (0,1,0,0), Y* = (0,0,1,0), Z" = (0,0,0,1). 

The expansion scalar 6 in the (0+l)D case reads 

e = A^d li u v = 

T 

The tensor a** 1 ' is also relatively simple when expressed in terms of the vectors X, Y, and Z 

a^(x) = W^U v) = oj_ (t) (X^X V + WW 1 ') + a^r) Z^ZY 
Here the functions a± and an are given by 


(26) 


(27) 


(28) 


a± = 


3r ’ 


3 T 


(29) 


The condition 2crj_ + try =0 reflects the fact that a ^ = 0. 


B. Phase-space distribution function for the (0+l)D case 

In the (0+l)D case, the anisotropy tensor in (8) reduces to the form 

= £_l(t) (Tw*' + y'W") 


-Ku { t ) z » z \ 


(30) 
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which is analogous to (28). Since is also traceless, the last equation reads 

= -^ (x»X v + WW") + £ Z^Z 1 ', 


(31) 


where we use the identification £ = . Thus, the distribution function in the local reference frame takes the form 

(32) 


f = N exp 


1 


—\\r 


+ (i + 0 - p! + (i + <l> + s) 


The form (32) suggests that it is useful to introduce the two new a variables defined as 

1 1 

a± = 


VI+ 0-e/2’ 

and to rewrite the LRF anisotropic distribution function as 


v/i + </> + £ 


/ = IV exp ( -- 




pi , a 

o' 0 

a, af. 


(33) 


(34) 


C. First-moment tensors 

Using the (0+l)D reduced anisotropic distribution function (34), the LRF energy density reads 


to 2 + ^=r + ] = ^3(a_L,aii,m)A 4 , (35) 

“x a \\ ] 

where rh = m/A and the function %3 has been defined in Ref. [40]. Because of the symmetry of the system the bulk 
viscosity II, being a scalar, must depend only on the proper time r, II = II(t). On the other hand, if we use the 
orthonormal basis presented above, we can write the shear tensor 7 jP" in the (0+l)D case as 

tt^" = ^tt s (t) ( X^X u + WW") - 7t s (t) Z^ZY (36) 

It turns out that it is completely described by the space basis {X,Y,Zj and the scalar function 7 t s (t). However, 
instead of using n s and n, it is sometimes convenient to use the longitudinal pressure Vl and the transverse pressure 
Vt, 

V L = Z-T-Z = V eq .+H-n s , V t = X-T-X = Y-T-Y = V eq .+n+K s , (37) 

since this results in relatively simple expressions, which do not depend on the effective temperature, namely 


£ = N d 3 p \/m 2 


p 2 exp 


A 


•Pi 


Vt 


N / d 3 p 




x/m^+p 2 6XP 1 A\ 


m 2 + P ± + P l 

«i a \\ 


~ /d 3 p 


p! 


A\ 


2 P 2 ± P\\ 

m 2 + Of- H-A 


= ‘H3L(ax,a||,m)A 4 , 
= H3T(o:± 1 ai\ 1 rh)X i . 


(38) 

(39) 


Here the integrals are taken in the local rest frame, and the functions 'Kzl and 77 st are also defined in Ref. [40]. 

One can recover the original degrees of freedom n and 7r s from the longitudinal pressure, the transverse pressure 
and the hydrostatic pressure in Eq. (4) using 


n = ^ (2 V T +V l )~ V eq ., ^s = \(Vt-V l ). (40) 

We note that only the integrals in Eqs. (35), (38), and (39) depend on the anisotropic distribution function. The rest 
of the arguments stem directly from the symmetry of the system, which is completely general in the (0+l)D set up. 
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D. Second- and zeroth-moment tensors 

Similarly to Eqs. (28) and (31), the tensor can be decomposed as 

Q pv = 0j_(t) ( 'x p X v + Y p yA + ©|| (t)Z p Z v , (41) 

and, consequently, we find 

0 <M = I ^ 0J _ _ (x p X v + Y p Y v ^j + | (©|| - 0 ± ) Z p Z v . (42) 

Using Eq. (34), the functions 0_l and ©y can be written as follows 

0_l = aa5_(47riV)A 5 TO 3 A'3(m), (43) 

©II = aa|(47T-ZV')A 5 m 3 A’3(m), (44) 

where a = a[]_ay. 

The last integral over the anisotropic distribution function which is necessary is the (non-equilibrium) number 
density 

n = an eq .( A), (45) 

with n eq (A) being the local equilibrium number density with the substitution T —> A. The equilibrium density reads 

7T eq . = (4 tt N)T 3 rhl q K 2 (rh eq .). (46) 

IV. REDUCED EVOLUTION EQUATIONS 
A. Energy conservation 

Since we are using the same form of the anisotropic distribution function as that used in Ref. [40], our equations 
expressing the energy-momentum conservation are the same as in [40]. Using Eq. (35), the Landau matching (10) can 
be written as 


7i 3 (a_L,ay,TO)A 4 = £ eq _(T,rh eq _). (47) 

The three-momentum conservation equations in (13) are automatically satisfied by symmetry. Indeed, using Eqs. (22) 
and (24) one finds that 

D = U-d = d T , Z-d=-d v , DU P = 0, (Z ■ d)Z p = -U p , d-Z = 0. (48) 

T T 

Using these relations we check that both the left- and right-hand sides of Eq. (13) vanish. In particular, for the 
right-hand side of Eq. (13) we obtain 


V #i 'P eq . ( T ) + V A1 n(r) + A p d a n pa = it a {r)Md a 


(x p X a + Y p Y a ^j - Z p Z° 


= —7t s (t) A p (Z-d)Z p + Z p d-Z 


= 0 . 


(49) 

(50) 


In a similar manner, from the energy conservation equation (12) one obtains 

jjg _ £ + P eq. + n — TTs _ £ + 'PL 


(51) 


Here we used the definition of the longitudinal pressure (37). Using Eqs. (35), (38), and (39) the last equation can be 
rewritten as 

(^'H 3 (a± , ay, m) — f2 m (aj_, ay, m)^ D In A + 2Q-r(a_L, ay, rh)D lnaj_ + Ui(a:j_, ay, hi) ^—b Z) In ay ^ =0. (52) 


For the explicit expressions of the functions fl m , Clr, and CIl see Ref. [40]. 








B. Shear viscosity equations 


The equations describing the shear pressure corrections are the same for the three different options for the bulk 
viscous pressure evolution equations, and they can be rewritten as 


K’ffDQ 010 + ^0^+2ejfV> A - 2Q^U} V)X = —— 0^0. 

o 


(53) 


' eq. 


Since the convective derivatives of X, Y, and Z either vanish or are proportional to U, the first term of the last 
equation reads 

A£A£D0 a/3 - * A^A a0 D@ a P = * (d©j_ - D0y) (x fl X u + WW") + (D0y - DQ ± ) Z»Z V . (54) 

The contraction of 0 with the expansion tensor a reads 

©i^ = ~\ (©KTI - ©m) {x»X v + w) - l ( ©||cr|| - ©idi) Z»Z\ (55) 

where we used the (0+l)D form of the expansion tensor (28). Since the vorticity w' 11 ' vanishes in the case of (0+l)D 
expansion, collecting all of the terms together in Eq. (53) one finds 


DQ ± - £>©„ + (© ± - ©||) ^ 2 (©xn - ©mi) + (©x - ©||) (x^X” + WW") + 


' eq. 


£>©x - £>©|| + -0 ^©j_ — ©||) — 2 ^©i<ti - ©|| CT li) + T — (©x — ©||) 


Z»Z V = 0. 


The equations for the shear pressure corrections are then straightforwardly fulfilled if 

DQ± — DO || + -6 ^0j_ — ©||) - 2 ^©i<Ti — ©||cr||) + —— ^0j_ — ©||) = 0- 


(56) 


(57) 


We note that the number of independent equations in the (0+l)D case is reduced to one. This is not surprising, 
since the number of equations matches the number of independent parameters defining the anisotropy tensor . To 
simplify the last expression we can use the explicit expressions for the four velocity gradients (29) 

( 0± _ 0|| ) + i (e ± - 0y) - *©y + ^ (©x - ©II) = 0. (58) 


D 


Using the expressions (43) and (44), and dividing (58) by the positive quantity a(4TrN)X 2 m 3 K 3 (m/X), Eq. (58) is put 
into the form 

K 2 ( to )' 


^4a j_ ~ 2a|) D In a± + — 3ay) D In ay + (a]_ — a|) 

(« 


5 + m 


K 3 ( to ) 


D In A 


(59) 


+ I aq — ax 


1 1 

- H- 

r r„ r 


2an 


= 0. 


C. Bulk viscosity equations 

1. Trace equation 

Using the decomposition (41) and (28) the trace of the second moment of the Boltzmann equation reads 


- D 


^2©x + ©||) — - 0 ^2©_l + 0||) +2 ^2©xcrj_ + ©|| <T ||) 


1 


3© eq . - (2©x + ©| 


(60) 


As was done for the shear pressure corrections equations, it is possible to use the Eq. (43), (44), (27), and (29) to 
further simplify the Eq. (60). Dividing by the positive function a(47rlV)A 2 m 3 K 3 (m/X) one obtains 

K 2 {m) 


— ^8a^ + 2a|) Din ax — ^2a^_ + 3ajj) Din ay — ^2a^ + ay) 

_( 9 „2 , A (l , J_\ S , 3 T 2 K 3 (m eq ) 

V ± Dvr r ea . j r r ea . aX 2 K 3 (m) 


5 + m 


A' 3 (m)_ 


D In A + 


(61) 
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2. U projection 

In the case where the U projection is applied we use the equation 


to A 2 (jti) < D In a 


Ki(rh) 

1 A 2 (to)_ 


I? In A 


1 

T 


1 


1 - 


+A 3 (to) ^ 8 a]^ + 2a|^ D lnaj_ — + 3a|^ D In ay — ^ 

-(2ai+t«|)(i+A-) 


T /v 2 (m eq .) 

aA A'2 (to) 


( 2 a | 4 - ai 


2 o 2 

**ll 

r 


3 T 2 /v 3 (m eq .)| _ 


a A 2 A 3 (to) 


= 0 . 


(62) 


5 + TO 


A 2 (to) 

A 3 (to)_ 


D In A+ 


5. Zeroth moment equation 

The zeroth moment of the Boltzmann equation was previously obtained in Ref. [40] using the same form (34) for 
the anisotropic distribution function, with the result being 


Aina + 


3 + TO 


Aj (to) 

A 2 (to) 


D In A H- 1 - 

T T eq . 


1 


T A 2 (TO eq .) 
aA A 2 (to) 


(63) 


V. NUMERICAL RESULTS 

In this section we present our numerical results. We first present comparisons of solution of the full non-perturbative 
equations of motion and then make comparisons using the linearized equations obtained in the near-equilibrium limit. 
In all figures, we take the initial proper time to be To = 0.25 fm/c and, when a fixed mass is required, we take 
to = 300 MeV. We evolve the equations subject to the assumption of fixed ratio of the viscosity, 77 , to equilibrium 
entropy density, S eq ,. The equilibrium entropy density S cq , at a given effective temperature T is obtained using the 
equilibrium relation 


Soq.(T) 


£eq .(T)+V eq .(T) 
T 


(64) 


while, for the shear viscosity, we use the formula [55, 56] 


V(T) = 


1 eq, 


(T)V eq (T) 


15 


eq. 


3 R 3 

TO 2 q A 2 


‘'eq. 


Ai _ Aip 
K-2 A 2 


(65) 


which is valid for a relativistic classical gas of massive particles. Here all the functions should be evaluated at m eq ., 
K n are modified Bessel functions, and 2Aj i = 7 r [1 — rh eq .Ko(rh eq .)L-i(m eq .) — rh eq .I\i(m eq .)Lo{rh eq .)] where Li is a 
modified Struve function [49]. From Eqs. (64) and (65) we find the relaxation time r eq . as a function of the effective 
temperature. The latter can be used directly in both hydrodynamic and kinetic calculations. In the remainder of this 
section, we will drop the ‘eq.’ subscript from the entropy with the understanding that the entropy appearing in the 
ratio r]/S is the equilibrium entropy density. 


A. Non-perturbative case 

In Fig. 1 we show our numerical results for the ratio of the longitudinal and transverse pressures, Vl/Vt , plotted 
as a function of the longitudinal proper time r. The behavior of this ratio is determined mainly by effects connected 
with the shear viscosity. The ratio determined from the various hydrodynamic calculations is normalized to the exact 
result obtained from kinetic theory . 1 The upper, middle, and lower panels of Fig. 1 are the results obtained with the 


1 For the detailed discussion of the method used for this purpose we refer the reader to [49]. 
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FIG. 1: (Color online) Ratio of the longitudinal and transverse pressures plotted as a function of the longitudinal proper time. 
The results of various hydrodynamic calculations are normalized to the exact result obtained from the kinetic theory. Details 
of the figure are explained in the text. 


initial anisotropy parameter £o = £(tq) = 0, 10, and 100, respectively. Similarly, the three columns are the results 
obtained using different values of the ratio of the viscosity to the entropy density; Atttj/S = 1 in the first column, 
47777/5 = 3 in the second column, and Anr/fS = 10 in the third column. 

The solid (brown), dashed-dotted (green), and dashed-double-dotted (violet) lines in Fig. 1 show the results obtained 
using leading-order anisotropic hydrodynamics with the use of Eqs. (61), (62), and (63), respectively. The long-dashed 
(red) line represents the results obtained with the recent second-order hydrodynamic formulation by Jaiswal [23], while 
the short-dashed (blue) lines show the results obtained within an earlier formulation of the leading order anisotropic 
hydrodynamics, designed for the (1+1)D massless case in [39] and generalized to the massive case by Nopoush et al. 
in Ref. [40]. The results presented in Fig. 1 indicate that, as long as the initial anisotropy is relatively small, £0 < 10, 
all versions of leading-order anisotropic hydrodynamics give better agreement with the exact result than second-order 
hydrodynamics equations of Ref. [23] . For larger initial anisotropies the versions of anisotropic hydrodynamics based 
on the trace or on the U projection of the second moment, see Eqs. (61) and (62), seem to be preferred. 2 


2 One should note here that, although the results obtained within the (1+1)D formulation obtained in Ref. [40] work better than all 
formulations obtained in [41] in the entire ( 77 , £ 0 ) parameter space, in their present form they cannot be used in the most realistic full 
































11 






r [fm/c] 


r [fm/c] 


r [fm/c] 


FIG. 2: (Color online) Proper time dependence of the bulk viscous pressure multiplied by the proper time. The hydrodynamic 
calculations are compared to the exact solution of the Boltzmann equation (thin solid black line). The remaining notation is 
the same as in Fig. 1. 


In Fig. 2 we show our numerical results for the proper time dependence of the bulk viscous pressure multiplied by 
the proper time. The notation is the same as in Fig. 1. In this case, the best description is given by the second-order 
hydrodynamics equations of Ref. [23], but anisotropic hydrodynamics using the zeroth moment, see Eq. (63), also gives 
satisfactory results. It is interesting to observe that very good agreement with the exact results is always obtained 
within the leading order (1+1)D formulation defined in Refs. [39, 40] (short-dashed blue lines). 


B. Close-to-equilibrium limit 

Despite the good agreement with the exact solutions seen in the last section, one may ask how much this agreement 
depends on the non-perturbative treatment of the large pressure anisotropy. It is possible that the close-to-equilibrium 
behavior is not well reproduced, but the non-linear treatment of the anisotropy compensates for it. In such a case, 
the set of equations would not be reliable for small pressure corrections. In heavy-ion collisions, however, the pressure 
corrections in the transverse plane are relatively small, therefore, the proposed equations should be reliable even for 


(3+l)D evolution. 
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FIG. 3: (Color online) Comparison of the transport coefficients. Left panel: Transport coefficients r 7!ltl <5^^, and A^-n from 
Ref. [23] (red lines) and [41] (black lines). Right panel: Transport coefficients An* (multiplied by 10) and <5nn from Ref. [23] 
(dashed lines) and from [41] (three different formulations shown as solid, dashed-dotted and dashed-double-dotted lines). All 
transport coefficients are plotted as functions of the ratio m/T. 




FIG. 4: (Color online) Left panel: Transport coefficients <f>e (solid black line) and fa (dashed red line) appearing in the evolution 
equation for the shear tensor in the linearised version (epe and fa are independent of the prescription chosen to describe the bulk 
pressure). Right panel: Transport coefficients fa (thick lines) and fa (thin lines) for three prescriptions of the bulk pressure 
evolution defined by Eqs. (61), (62) and (63). 


small pressure corrections. 

It has been shown in Ref. [41] that in the close-to-equilibrium limit all prescriptions for the (3+l)D anisotropic 
hydrodynamics dynamical equations give a structure that is of the same form as obtained using second order viscous 
hydrodynamics [24]. One can easily show that, for the Bjorken scenario considered herein, Eqs. (42) and (43) from 
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FIG. 5: (Color online) Same as Fig. 1 but for the close-to-equilibrium limit. 


Ref. [41] reduce to 


DU 


n C A n 

— =-oim — 

Tfr TnT T 


H-= -- 

Ttt OT-jrT 


1 

l 


r r n r n 2 
\ 7T S 2 n , fa ^ n fa 7T 2 

^7T7T ) I 0 A^n I TTfi-l-J- n • 

T 6 T T^r TV 2 


( 66 ) 

(67) 


They have to be considered together with Eq. (51). The transport coefficients do not depend on the anisotropy or the 
symmetry of the expansion. We can therefore directly compare the transport coefficients derived in [41] with other 
results presented in the literature. This is done in Figs. 3 and 4. 

In the left panel of Fig. 3 the coefficients T n7r ,8 nn and A^n from Ref. [23] and [41] are plotted. They appear in 
the linearised equation for the shear tensor. We observe signficant differences between the results derived in [23] and 
[41], in particular, for the coefficient A^n, which diverges in the high-temperature limit of the linearised version of 
anisotropic hydrodynamics. In the right panel of Fig. 3 the transport coefficients An^ and <5nn from Ref. [23] and 
from [41] are plotted. They appear correspondingly in the linearised version of the equation for the bulk pressure. 
Again, one observes noticeable differences between the results derived in [23] and [41]. 

In the left (right) panel of Fig. 4 we present the kinetic coefficients fa and fa (fa and fa) used in the linearised 
equation for (II). All these results follow from the linearised version of anisotropic hydrodynamics [41]. In the 
analogous second-order hydrodynamics formulations all the (f> coefficients vanish for the collision term treated in the 
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FIG. 6: (Color online) Same as Fig. 2 but for the close-to-equilibrium limit. 


relaxation time approximation considered in [23]. 

Interestingly, despite of the differences between the kinetic coefficients obtained in [41] and [23], the final results 
for the pressure anisotropy and the bulk pressure are quite similar, see Figs. 5 and 6 , particularly if one uses the 
zeroth moment in the bulk pressure evolution equation. This behaviour may be connected with the values of the 
kinetic coefficients </>i and </>3 which have noticeably smaller values for the zeroth moment case. It is also interesting 
to compare Fig. 1 with Fig. 5, and Fig. 2 with Fig. 6 . This comparison indicates that it it is generally necessary to 
take into account the full equations of anisotropic hydrodynamics in order to have the best agreement with the exact 
solution of the kinetic equation. 


VI. SUMMARY AND CONCLUSIONS 

Our results show that the use of the recently formulated leading-order anisotropic hydrodynamics is comparable 
with the present best prescriptions of the standard second-order viscous hydrodynamics. In specific applications, we 
find that the behavior of the bulk viscous pressure is better described if the zeroth moment of the kinetic equation is 
used, while the ratio of the longitudinal to transverse pressure is better reproduced if one uses the trace of the second 
moment or the projection of the second moment on the flow vector. Our results, taken together, suggest that the 
overall best option for the (3+l)D anisotropic hydrodynamics formulated in Ref. [41] is to use the zeroth moment 
of the Boltzmann equation to derive the bulk viscous dynamics. Looking at both the close-to-equilibrium and fully 
non-linear results we observe that the use of the non-linear version generally improves the agreement with the exact 
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solutions. This also suggests that the equations derived in Ref. [41] may be successfully used as the first term in the 
expansion of (3+l)D anisotropic hydrodynamics. 
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